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Tests of CPT and Lorentz symmetry using neutral-meson oscillations are studied within a for- 
malism that allows for indirect CPT and T violation of arbitrary size and is independent of phase 
conventions. The analysis is particularly appropriate for studies of CPT and T violation in os- 
cillations of the heavy neutral mesons D, Bd, and B s . The general Lorentz- and CPT-breaking 
standard-model extension is used to derive an expression for the parameter for CPT violation. 
It varies in a prescribed way with the magnitude and orientation of the meson momentum and 
consequently also with sidereal time. Decay probabilities are presented for both uncorrelated and 
correlated mesons, and some implications for experiments are discussed. 
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I. INTRODUCTION 

The original discovery of CP violation in the neutral- 
kaon system jjj has led to numerous theoretical and 
experimental studies of discrete symmetries in neutral- 
meson oscillations 0j. Much of the effort has been fo- 
cused on the K system, but the advent of high-statistics 
experiments involving the heavy neutral mesons, in par- 
ticular the Bd mesons ||, has opened the door for a 
broader class of investigations. 

In a neutral-meson system, the violation of CP sym- 
metry includes the possibility of CPT violation H,|). 
For the K system, CPT violation in oscillations can be 
parametrized by a complex quantity 6k that is known to 
be small or zero Q. Under the ad hoc assumption that 
6k is a constant complex number, experiments have es- 
tablished that its real and imaginary parts are no greater 
than about 10" 4 0§. 

The assumption of constant nonzero 6k is known to 
fail in conventional quantum field theory. The nature 
of 6k is determined by the properties of the theory un- 
der Lorentz transformations. For any realistic Lorentz- 
invariant quantum field theory such as the standard 
model, the CPT theorem shows that 6k must be zero [|J . 
If instead Lorentz violation is allowed, then using an ex- 
plicit and general standard-model extension || to calcu- 
late 6k reveals that it varies with the meson 4-momentum 
^0| , p| . This variation has recently been exploited by the 
KTeV Collaboration in placing a qualitatively new bound 
on CPT violation in the neutral- if system 

For systems involving the heavy mesons D, Bd, B s sev- 
eral CPT tests have been proposed ]T^-p^|, and bounds 
have been obtained in some recent experiments with the 
Bd system All these results rely on the assumption 
of a nonzero constant complex parameter for CPT viola- 
tion. However, as in the K system, this assumption fails 
in realistic quantum field theories: either the parame- 
ter vanishes by the CPT theorem, or it depends on the 
4-momentum of the meson. 

The present work provides a general treatment of CPT 



violation in neutral-meson oscillations in the context of 
quantum field theory allowing for Lorentz violation. A 
convenient formalism is adopted that is independent of 
phase conventions and allows for CPT and T violation 
of arbitrary size in any neutral-meson system. The com- 
plex parameter for CPT violation is calculated in the 
general Lorentz-violating standard-model extension, re- 
vealing a well-defined variation with the magnitude and 
orientation of the meson momentum and a correspond- 
ing variation with sidereal time. Some experimentally 
relevant decay probabilities and asymmetries are derived 
for both uncorrelated and correlated mesons. The results 
obtained here complement the analyses in earlier works, 
which described some essential physics @ and obtained 
expressions valid for small CPT violation in the K, D, 
Bd, and B s systems 

Section |0] provides background information and fixes 
some notational conventions. A suitable parametrization 
of the effective hamiltonian for the time evolution of a 
neutral-meson state with CPT and T violation of arbi- 
trary size is presented in section III. The calculation 
of the complex parameter for CPT violation is given in 
section [IVL Implications for experiment are considered 
in section |Vj. The appendix contains a brief description 
of other formalisms adopted in the literature. Through- 
out this work, a strong-interaction eigenstate is denoted 
generically by P°, where P° is one of K°, D°, B° d , B° s , 
and the corresponding opposite-flavor antiparticle is de- 
noted P°. 



II. BASICS 



An arbitrary neutral-meson state is a linear combina- 
tion of the Schrodinger wave functions for the meson P° 
and its antimeson P°. This combination can be repre- 
sented as a two-component object ^(t), with time evolu- 
tion governed by a 2 x 2 effective hamiltonian A according 
to the Schrodinger-type equation M 
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idtV = A*. (1) 

Throughout this paper, subscripts P are understood on 
'J, on the components of the effective hamiltonian A, and 
on related quantities such as meson masses and lifetimes. 

The physical propagating states are the eigenstates of 
A, analogous to the normal modes of a classical two- 
dimensional oscillator jTjJ . In this work, these states are 
generically denoted as \P a ) and \P b ). They evolve in time 
as 

\P a {t)) = exp{-iX a t)\P a ), 
\P b {t))=exp{-iX b t)\P b ). (2) 

The complex parameters A a , X b are the eigenvalues of A. 
They can be decomposed as 

A a = m a - i«7 a , X b = m b -\i^ b , (3) 

where m a , m b are the propagating masses and j a , j b are 
the associated decay rates. For the K system, contact 
with the standard notation can be made via the identifi- 
cation m a = ms, m b = m L , j a = 75, "f b — j L . For the 
D system, there is no well established convention and I 
use the notation in Eq. (||). For the and B s systems, 
the relation to the standard notation can be taken as 
m a = m L , m b = m H , j a = T L , j b = T H . 

For calculational purposes, it is useful to introduce a 
separate notation for the sums and differences of these 
parameters: 

A = A Q + X b = m - |vy, 
AA = A a - A 6 = -Am - ±iA 7 , (4) 

where m — m a + m b , Am = m b — m a , 7 = 7„ + 75, 
A7 = 7a — 76- Note in particular the choice of sign in the 
definition of A7, which coincides with that in the K sys- 
tem but is the negative of the quantity Ar often adopted 
in the P>d system. The reader can therefore make direct 
contact with results in the latter convention by identify- 
ing A7 = — Ar in any equation in this work. 

The off-diagonal components of A control the flavor 
oscillations between P° and P°. Indirect CPT violation 
occurs if and only if the difference of diagonal elements of 
A is nonzero, An — A22 7^ 0. Indirect T violation occurs 
if and only if the magnitude of the ratio of off-diagonal 
components of A differs from 1, IA21/A12I 7^ 1. 

A priori, the effective hamiltonian A can be 
parametrized by eight independent real quantities. Four 
of these can be specified in terms of the masses and decay 
rates, two describe CPT violation, and one describes T 
violation. The remaining parameter, determined by the 
relative phase between the off-diagonal components of A, 
is physically irrelevant. It can be dialed at will by rotat- 
ing the phases of the P° and P° wave functions by equal 
and opposite amounts. The freedom to perform such ro- 
tations exists because the wave functions are eigenstates 



of the strong interactions, which preserve strangeness, 
charm, and beauty. Under a rotation of this type involv- 
ing a phase factor of exp(ix) for the P° wave function, the 
off-diagonal elements of A are multiplied by equal and op- 
posite phases, becoming exp(2«x)Ai2 and exp(— 2ix)^2i- 

III. FORMALISM 

Since relatively little experimental information is avail- 
able about CPT and T violation in the heavy neutral- 
meson systems, a general parametrization of A is appro- 
priate. It is desirable to have a parametrization that 
is model independent, valid for arbitrary size CPT and 
T violation, independent of phase conventions, and ex- 
pressed in terms of mass and decay rates insofar as pos- 
sible. A parametrization of this type was originally in- 
troduced by Lavoura in the context of the kaon system 
PUP^f . For simplicity, it is also attractive to arrange 
matters so that the quantities controlling T and CPT vi- 
olation are denoted by single symbols that are distinct 
from other frequently used notation. In this section, a 
parametrization convenient to the four meson systems 
and satisfying all the above criteria is presented and re- 
lated to formalisms often used in the literature. 

For a complex 2x2 matrix, it is possible to write the 
two diagonal elements as the sum and difference of two 
complex numbers. It is also possible to write the off- 
diagonal elements as the product and ratio of two com- 
plex numbers. Using these two facts, which ultimately 
permit the clean representation of T- and CPT- violating 
quantities, a general expression for A can be taken as: 

(u+i vw- 1 \ 

A = iAA , (5) 

\ VW U - £ / 

where the parameters UVW^ are complex. The factor 
A A/2 has been extracted from A to make these parame- 
ters dimensionless and to avoid factors of 2 in expressions 
below. 

The requirements that the trace of the matrix is tr A = 
A and that the determinant is det A = X a X b impose the 
identifications 

U=X/AX, V=^l-e (6) 

on the complex parameters U and V. The free pa- 
rameters in Eq. (||) are therefore W and £. These can 
be regarded as four independent real quantities: W = 
wexp(iw), £ = Re£ + ilm£. One of these four real num- 
bers, the argument lu of W, is arbitrary and physically 
irrelevant. It changes under the phase redefinitions dis- 
cussed at the end of the previous section. The other three 
are physical. The modulus w of W controls T violation, 
with w = 1 if and only if T is preserved p0|. The two 
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remaining real numbers, Re£ and Im£, control CPT vio- 
lation and both are zero if and only if CPT is preserved. 
The quantities w and £ can be expressed in terms of the 



components of A as 21 



= VlWAul, £ = AA/AA, 



(7) 



where AA = An — A 22 . 

In this wt; formalism, the three parameters for CP vio- 
lation w, Re£, Im£ are dimensionless and independent of 
phase conventions. They are phenomenologically intro- 
duced and therefore are independent of specific models. 
However, this does not imply that they are necessarily 
constant numbers. Indeed, the assumption of constancy 
for £ frequently made in the literature is a special choice 
that strongly restricts the generality of the parametriza- 
tion and which according to the CPT theorem is inconsis- 
tent with the fundamental structure of Lorentz-invariant 
quantum field theory. In fact, if the requirement of exact 
Lorentz symmetry is relaxed, then £ cannot be a constant 
quantity within the framework of quantum field theory 
but instead must vary with the momentum of the meson. 
Since CPT violation is a profound effect, it is unsurpris- 
ing that the parameter £ has features different from w. 



The choice of the notation £ (rather than, say, X) in Eq. 
(||) has been made to emphasize this crucial fact. 

The physical states with definite mass and lifetimes are 
the eigenstates of A. In the u>£ formalism, they take the 
form 



\Pa) = M a (\P°) + A\P°)), 

\P b )=M b (\P°)+B\n), 



where 



A = (1 - i)W/V, B = —(1 + CjW/V. 



(8) 



(9) 



The normalizations Af a , ■N'b m Eq. (j|) can be chosen as 
desired. For unit-normalized states, the normalizations 



arc 



(ir, a )/VT+\Aj 



Af a = exp 



N b = exp(i»fr)/vi + \B\ 



(10) 



where r\ a and rj b are phases that can be chosen freely. 
For the analysis of physical observables in the following 
sections, the values of these phases are irrelevant p2|. 



Formalism 


Parameters depend 
on phase convention? 


A, AA 
given as 


CPT parameter 
(complex) 


T parameter 
(real) 


wt; 


No 


A, AA 




w 


MY 


Yes (M 12 , r 12 ) 


See Eq. @ 


(M n ~ M 22 ) 
-Wii-T 22 ) 


\M 12 -ir 12 /2\ 




Yes (E u E 2 ) 


-2W, 2^jE\ + E% + El 


E3 


i(EiE* - E{E 2 ) 


DE9<j) 


Yes (0) 


-2iD, 2E 


cos 9 


|exp(i0)| 


pqrs 


Yes (p, q, r, s) 


A, AA 


(ps - qr) 


\pr/qs\ 


eS 


Yes (e, 6) 


A, AA 


S 


Re e, if CJP small 



Table 1: Comparison of formalisms for neutral- meson mixing. 



Some insight into the advantages of the w£ formal- 
ism can be obtained by comparing it to alternative for- 
malisms available in the literature. The appendix sum- 
marizes some of the more popular ones, and Table 1 pro- 
vides a comparative synopsis of their features. The first 
column identifies the formalism through the standard no- 
tation for its parameters. The second column indicates 
the phase-convention dependence of its parameters. The 
third column lists the connection between the physical 
quantities A, AA and their expression in the given for- 
malism. The fourth column specifies the complex combi- 
nation of parameters that governs CPT violation in the 
specified formalism, while the last column gives the real 
number controlling T violation. Note that the final entry 
on the last line holds only for small CPT and T violation 



and assumes a phase convention with Im e = 0. 

Exact relationships exist between the wt; formalism 
and the other formalisms listed in Table 1, but they can 
be involved and may change with the choice of phase con- 
ventions. Expressing the complex parameter £ for CPT 
violation in the other parametrizations gives 



£ = %[(Mu - M 22 ) - |i(T u - T 22 )] 

x[(M 12 -|ir 12 )(M* 2 -|irj 2 ) 

1 [(A^ii-^2 2 )-^(rii-r 22 )] 2 ] 



4 LV— J H 

E 3 



-1/2 



= cos 9 



E 2 2 



El 



3 



(ps - qr) 
(ps + qr) 
28. 



(11) 



The last line is valid only for small 8 and e and only in a 
special phase convention, but shows that £ can be identi- 
fied with 28 for an appropriate choice of phase convention 
in the K system. In any case, for the D, Bd, and B s sys- 
tems, £ appears simpler to use than 8 or any of the other 
parametrizations . 

A similar exercise for the real parameter w for T vio- 
lation yields 
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= \(M* 2 -^r* 12 )/(M 12 -^r 12 )\^ 

= \(E 1 +iE 2 )/(E 1 -iE 2 ) 
= |exp(#)| 

= V\qs/pr\ 

w 1 - 2Ree. 



(12) 



The last line is again valid only for small 8 and e and 
only in a special phase convention. 

The above equations reveal that the wf; formalism is 
most closely related to the DE6<p formalism, but offers 
a more direct link to A, AA, an abbreviated notation 
for CPT violation, and a single symbol for the phase- 
independent physical parameter for T violation. On the 
more practical side, the use of £ also avoids confusion 
with the standard use of the track orientation angles 9, 
<j> for the meson in the detector, which is a useful asset in 
the presence of orientation-dependent CPT-violating ef- 
fects. Overall, advantages of the wt; formalism include its 
model independence, its use of mass and decay rates as 
physical parameters, its validity for arbitrary-size CPT 
and T violation, and its independence of phase conven- 
tions. In the present work, use of the u>£ formalism sim- 
plifies the results of the study of CPT violation. 



IV. THEORY FOR CPT VIOLATION 

The CPT theorem guarantees CPT invariance of 
Lorentz-symmetric quantum field theories, including the 
usual standard model of particle physics. To construct a 
description of CPT violation viable at the level of quan- 
tum field theory, it is therefore of interest to consider 
the possibility of small violations of Lorentz invariance. 
A general standard-model extension allowing for Lorentz 
and CPT violation is known M. It could emerge, for ex- 
ample, as the low-energy limit of a fundamental theory 
at the Planck scale p3[ . This standard-model extension 
provides a quantitative microscopic theory for Lorentz 
and CPT violation that is applicable to a wide class of 
experiments in addition to the studies of neutral-meson 
oscillations considered in the present work. Among these 
are, for example, comparative tests of QED in Penning 



tra ps |24| - |27[ , spectroscopy of hydrogen and antihydrogen 
p8| , [29|] , measurements of muon properties ||(]j3l| , clock- 
comparison experiments [|32"|-[35|| , observations of the be- 
havior of a spin-polarized torsion pendulum J36|,^7t, mea- 



surements of cosmological birefringence Hggggg, and 
observations of the baryon asymmetry |4l|. However, 
none of these tests are sensitive to the sector of the 
standard-model extension involved in the experiments 
with neutral-meson oscillations, essentially because the 
latter are flavor changing JlO[ . 

Using the general standard-model extension, a per- 
turbative calculation can be performed to obtain the 
leading-order CPT-violating contributions to A. These 
emerge as the expectation values of interaction terms 
in the standard-model hamiltonian ||L_ The CPT- 
unperturbed wave functions \P°) and \P°) are the ap- 
propriate states for constructing the expectation values. 
The hermiticity of the perturbation hamiltonian ensures 
reality of the dominant contributions to the difference 
AA = An — A22 of the diagonal terms of A and therefore 
constrains the form of A. It can be shown that [[To) 



AA w ^Afl, 



(13) 



where (3^ = 7(1, (3) is the four- velocity of the meson state 
in the observer frame. The effect of Lorentz and CPT vi- 
olation in the standard-model extension appears in Eq. 
( |l3| ) via the factor Aa^ — r qi a* 1 — r q2 ajf , where a -? 1 , a* 2 
are CPT- and Lorentz-violating coupling coefficients for 
the two valence quarks in the P° meson, and where r qi 
and r q2 are quantities resulting from quark-binding and 
normalization effects . The coefficients a* 1 , a* 8 for 
Lorentz and CPT violation have mass dimension one and 
emerge from terms in the lagrangian for the standard- 
model extension of the form — a^<77 M q, where q specifies 
the quark flavor. 

The 4- velocity and hence 4-momentum dependence in 
Eq. (|l^) confirms the failure of the usual assumption of 
a constant parameter for CPT violation. This depen- 
dence has substantial implications for experiments, since 
CPT observables will typically vary with the momentum 
magnitude and orientation of the mesons. As a result, 
the CPT reach of an experiment is affected by the meson 
momentum spectrum and angular distribution |ic|]ll|| . 

A significant consequence of the 4-momentum depen- 
dence arises from the rotation of the Earth relative to the 
constant vector Aa. This leads to sidereal variations in 
some observables juJIyJ. The point is that the analysis 
leading to Eq. ( fl3| ) is performed in the laboratory frame, 
which rotates with the Earth. The resulting sidereal time 
dependence can be exhibited explicitly by converting the 
expression for AA to a nonrotating frame. 

Denote the spatial basis in the nonrotating frame by 
(X . Y . Z) and that in the laboratory frame by (x,y, z). 
Following Ref. p3], define the nonrotating- frame basis 
(X,Y,Z) to be compatible with celestial equatorial co- 
ordinates E3j with Z aligned along the Earth's rotation 
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axis. The z axis in the laboratory frame can be chosen for 
maximal convenience. For collimated mesons, it may be 
useful to take it as the beam direction. In a collider, the 
direction of the colliding beams could be adopted. For a 
nonzero signal involving sidereal variations, cosx = z- Z 
is nonzero, and z precesses about Z with the Earth's 
sidereal frequency f2. A complete map between the two 
bases is given by Eq. (16) of Ref. |Q. For convenience 
in what follows, take 9 and <f> to be conventional polar 
coordinates defined about the z axis in the laboratory 
frame. If the z axis is chosen along the axis of a detector, 
then 9, cj> are the usual detector polar coordinates. 

Any coefficient a for Lorentz violation with laboratory- 
frame components (a 1 , a 2 , a 3 ) has nonrotating- frame 
components (a x ,a Y ,a z ) given by Eq. (12) of Ref. Q |. 
This relation determines the sidereal variation of Aa and, 
using Eq. (jl^), of A A. The complete momentum and 
sidereal-time dependence of the parameter £ for CPT vi- 
olation in any of the P systems can then be obtained. 
Noting that the laboratory-frame 3-velocity of a P me- 
son has the form /3 = f3 (sin 9 cos 4>,sm 9 sin <fi, cos 9), and 
the momentum magnitude is p = \p\ = Pmp r y(p), where 
7(p) = y/1 + p 2 /nip as usual, the expression for £ is 
found to be 



AA 



{ Aao + (3Aaz (cos 9 cos \ — sin 9 cos <j) sin x) 

+/3 [Aay (cos 9 sin x + sin 9 cos 4> cos x) 

—Aax sin 9 sin </)] sin Qt 
+(3 [Aax (cos 9 sin x + sin 9 cos (j) cos x) 

+ Aay sin 9 sin <p\ cos ili } , (14) 



where i denotes the sidereal time. 

In deriving Eq. (|l4|), only leading-order terms in a M 
have been kept but no other assumption about the size 
of £ has been made. The result (O) is therefore a gen- 
eralization of Eq. (13) in Ref. |11[|, which was obtained 
for the K system under the assumption of small 5k ■ In 
particular, Eq. ([l4]) holds for the heavy-meson systems 
where the possibility of large |£| > 1 remains experimen- 
tally admissible at present. 

Note that the expressions ( |l3"| ) and (|l4j) explicitly show 
that the real and imaginary parts of £ are connected 
through the mass and lifetime differences of the two phys- 



ical eigenstates P a , Pb |13|. The relationship is 



Re£ = -2AmIm£/A7. 



(15) 



However, in the interest of generality this result is used 
only sparingly in this work. 



V. EXPERIMENT 

To illustrate some implications of the result ( fl4| ) , this 
section derives some experimentally relevant decay am- 
plitudes, probabilities, and asymmetries. For simplicity, 
attention is restricted to the case of semileptonic decays 
into a final state / or its conjugate state /. Although 
studying these decays suffices for present purposes, other 
decays are also likely to be relevant in practice, and it 
would be of interest to perform a more complete study. 
Another simplification adopted here is the neglect of any 
violations of the AQ = AS, AQ = AC, or AQ = AB 
rules. A careful consideration of these and other more 
mundane complications would certainly be important in 
a definitive experimental analysis [|l9|] . However, since 
there is no reason to expect such complications to exhibit 
observable momentum or sidereal-time dependences, the 
extraction of a compelling positive signal for CPT viola- 
tion should be feasible. 

Under these assumptions, the basic transition ampli- 
tudes for semileptonic decays can be taken as 



(f\T\P°)=F, (f\T\P°) = 0, 
(f\T\P°) = F, (f\T\P°) = 0. 



(16) 



Note that this parametrization allows for direct CPT vi- 
olation, which is proportional to the difference F* — F, 
as well as direct T violation. 

To determine the time-dependent decay amplitudes 
and probabilities, it is useful to obtain an explicit ex- 
pression for the time evolution of the neutral-P states. 
The wave functions \P°) and \P°) can be constructed in 
terms of \P a ) and |P&), and their evolution with the me- 
son proper time t can then be incorporated via Eq. (|^). 
This gives 



P°(t,i,P) 



c + se, svw 
svw^ 1 C-S£ 



The functions C and S depend on the meson proper time 
t and are given by 

C = cos(iAAi) exp(-iiAt) 



= w 



-i\ a t 



-iXht 



S = -i sin(i AXt) exp(-^iXt) 



|(e 



.2 
-i\ a t 



-IXht 



)■ 



(18) 



In addition to the proper-time dependence in S and C, 
Eq. (^) also contains sidereal time and momentum de- 
pendence from £(t,p). Since the meson decays occur 
quickly on the scale of sidereal time, it is an excellent 
approximation to treat sidereal time t as a parameter 
independent of the meson proper time t. It is therefore 
appropriate to take £ as independent of t but varying with 
t. This approximation is implemented in what follows. 
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A. Uncorrelated Mesons 

For the case of uncorrelated meson decays, the time- 
dependent decay probabilities can be obtained by com- 
bining Eqs. (Til) and (fTS). This gives 



p f (t,i,p)^\(f\T\p(t),i,p)\ 2 

x [(! + |C| 2 )coshA7t/2+ (1 - |£| 2 )cosAmi 
- 2Rc£sinhA7</2- 21m £ sin Amt] , 
Pj(t, t,p) = \(f\T\P(t, t,p))\ 2 =Pf(Z^ -f, F - F), 

Pj{t,t,P) = \{!\T\p{t,i,P))\ 2 

= \\F\ 2 w 2 \l - £ 2 | e-^/ 2 (cosh A 7 i/2 - cos Ami), 

l/tu,F->F), 
(19) 



p f {t,t,P) = \(f\T\p(t,i,m 2 = p 7 (w 



where the dependence on sidereal time t and momentum 
p is inherited from that of £ in Eq. (|l4|). Inspection of 
these equations reveals that nonzero indirect CPT viola- 
tion changes the shape of the first two probabilities, while 
both CPT and T violation merely scale the latter two. I 
emphasize that these expressions are valid for CPT and 
T violation of arbitrary size. They are also manifestly 
independent of the choice of phase convention H] . 

To extract the CPT and T violation from the decay 
probabilities (|l9|), it is useful to construct appropriate 
asymmetries. For the case of T violation, the dependence 
on sidereal time and meson momentum has relatively lit- 
tle effect. The last two probabilities in Eq. (|l9|) have the 
same CPT but different T dependences, and their differ- 
ence divided by their sum is sensitive to the parameter 
w for T violation but independent of the parameter £ for 
CPT violation and hence independent of sidereal time 
and meson momentum. In contrast, for the case of CPT 
violation the situation is more involved and several new 
features appear. 

As a simple example illustrating some of the effects, 
consider the case where F* = F, i.e., neglible direct 
CPT violation. The usual procedure is to assume con- 
stant nonzero £ (which is inconsistent with quantum field 
theory, as discussed above) and define an asymmetry 
-<4cpt(£) for CPT violation as 



P 7 (t)-P f (t) 
CPTU Pj(t) + P f (t) 



(20) 



The comparable definition in the present context is still 
useful but results in an asymmetry depending also on 
sidereal time and meson momentum: 



AcPT{t,i,p) 



Pj(tJ,p) - P f (t,t,p) 



P 7 (t,t,p) + P f (t,t,p) 



2Ref sinh A 7 i/2 + 2Imf sin Amt 
(1 + |£| 2 )coshA 7 i/2+ (1 - |f| 2 )cosAmt' 



(21) 



where the t, p dependence of f is understood. 

In practice, the efficient practical application of this 
and related asymmetries depends on the nature of the 
experiment. Appropriate averaging over one of more of 
the variables t, t, p, 9, 4> either before or after construct- 
ing the asymmetry ( |2l| ) can aid the clean extraction of 
bounds on Aa^. For instance, under certain circum- 
stances it may be useful to sum the data over <f> and 
use an asymmetry like Eq. (21) but defined with the (\>- 
average of Eq. (|l^). The form of Eq. (^J) shows that 
binning the data in t typically provides information on 
Aax and Aay, while binning in permits the separa- 
tion of the spatial and timelike components of Aa M . The 
p dependence can also be useful |]l0| , |lT| . 

As a specific example, already used in the K system 
PUP^f , suppose the mesons involved are highly colli- 
matcd in the laboratory frame. Then, the 3-velocity can 
be written (3 — (0,0, (3) and the expression ( |l4| ) for £ 
simplifies to 



+f3 sin x( Aay sin fit + Aa x cos . (22) 



Binning in t therefore provides sensitivity to the equato- 
rial components Aax , Aay , while averaging over t elim- 
inates them altogether. Indeed, a conventional measure- 
ment that ignores the dependence on sidereal time and 
meson momentum is typically sensitive only to the aver- 
age magnitude 



|f | =7|Aa + /3Aa z cosx|/|AA| 



(23) 



where (3 and 7 are averages weighted over the meson- 
momentum spectrum. This shows explicitly that previ- 
ous analyses performed under the assumption of constant 
CPT parameter produce results dependent on the type 
of experiment. 

If CPT violation is small so £ < 1 , the asymmetry d2l| ) 
takes the form 



AcvT{t,i,p) 



2Ref sinh A 7 i/2 + 2Imf sin Amt 
coshA7t/2 + cosAmi 



(24) 



A further assumption that could be countenanced in- 
volves the approximation of small Ajt/2, i.e., t < 2/A 7 . 
This gives 



A C PT{t,i,p) 



Re £ Ajt + 21m £ sin Amt 
1 + cos Amt 



(25) 



It is tempting also to neglect as small the term involv- 
ing Re£, but this is potentially invalid because Re£ ex 
Im£/A 7 according to Eq. (|l5|). Imposing the prediction 
(|H)|) instead gives 



G 



■Acpt (*,*>£) 



2Im£(sin Ami — Ami) 
1 + cos Ami 



(26) 



The extraction of complete information about Aa M re- 
quires clean CPT tests involving asymmetries such as Eq. 
( pl| ) that are independent of the parameter w for T viola- 
tion. However, the dependence on sidereal time of certain 
CPT-violating effects offers the possibility of extracting 
clean CPT bounds on spatial components of Aa^ even 
using observables that mix T and CPT effects jnj. An 
example is provided by the standard rate asymmetry 5i 
for Kl semileptonic decays (43): 



T(K L -> Z+7T-I/) - T(K L -> r-K+V) 
T(K L 1+tx-v) + T(K L -> l-Ti+v) ' 



(27) 



which under the assumption of constant nonzero param- 
eter for CPT violation (inconsistent with quantum field 
theory, as noted above) is determined by a combination of 
T and CPT effects that cannot be disentangled without 
further information. In the formalism, the asymme- 
try ( p7f ) and its generalization to arbitrary P& is found to 
be 



T(P b -> f) - T(P b -> /) 

r(p 6 -> /) + r(p -> /) 
|i-eV|i + £l 2 ^ 2 

\i - e\ + \i + z\ 2 w 2 

(l-ti;)- Re 



(28) 



where the last line assumes w ~ 1, £ *C 1, i.e., small T 
and CPT violation. Binning in sidereal time or momen- 
tum can therefore under suitable circumstances bound 
the spatial components of Aa^ independently of T vio- 
lation, even for observables involving both T and CPT 
violation. 



f 2 at proper time t 2 . As before, in tracking the sidereal- 
time dependence, it is an excellent approximation to re- 
gard the time interval between quarkonium production 
and detection of the decay products as negligible on the 
scale of the Earth's rotation period, so in what follows 
the creation of the state \i) and its evolution through the 
double decay process are taken to occur at fixed sidereal 
time t. 

The probability amplitude Af x f 2 for the double decay 
can be regarded as a function of the decay times t\, t 2 , 
of the sidereal time i, and of the two meson momenta pi, 
p 2 . It is given by 

A hh = Af 1 f 2 (t 1 ,t 2 ,t,p 1 ,p 2 ) = (hf2\T\i) 

= ^[(f 1 \T\P°(t u ij 1 ))(f 2 \T\P^(t 2 J,p 2 )) 

- (h\T\pZ(ti,i, p,)) (f 2 \T\P°(t 2 , t, p 2 ))] . (30) 

The time evolutions of \P°(t,t,pj) and \P°(t,t,p)) are 
determined by Eq. (|l7|). In substituting these expres- 



sions into the decay amplitude (30), care is required to 
keep separate track of the CPT-violating parameters £i 
and £2 for each meson, since they depend on the meson 
3-momenta and therefore typically differ in accordance 
with Eq. @. 

It is convenient and feasible to write a single expres- 
sion that holds for all double decay modes, including the 
various double-semileptonic combinations. For a = 1,2, 
define 



(fa\T\P°) - F a 



(fa\T\P°)=F a 



(31) 



and let C a = C(t a ), S a = S(t a ). Then, the probability 
amplitude is found to be 



B. Correlated Mesons 

Another situation of experimental importance is the 
case of correlated meson pairs, resulting from quarko- 
nium production and decay. The normalized initial quan- 
tum state ensuing immediately after the strong decay of 
the quarkonium can be written as 

\i) = -L(|P°(+))|P0(-)) - |P°(-))|po (+ ))), (29) 

where (+) indicates the meson travels in a specified di- 
rection in the quarkonium rest frame while (— ) indicates 
it travels in the opposite direction. Note that this initial 
state is independent of the choice of phase convention. 

Let the meson moving in the (+) direction have 3- 
momentum p\ in the laboratory frame and decay into a 
final state f± at proper time t\. Similarly, let the other 
meson have 3-momentum p 2 and decay into a final state 



+ (F 1 F 2 - F 2 F 1 )(C l C 2 - (£16 + V 1 V 2 )S 1 S 2 ) 

+ (F1F2VF- 1 - PiP 2 W0(y 2 CiS 2 - ViSiCh) 

+ (P1P2VK- 1 +F 1 F 2 W){^V 2 -bVi)SiS 2 ], (32) 

where the dependence on t and pi, p 2 is understood. The 
quantities V\, V 2 are defined in terms of £1, £2 by Eq. (||), 
while W — wexp(iuj) as before. 

Next, consider the special case of double-semileptonic 
decays and adopt the notation of Eq. ( pif ) . It is useful to 
introduce the definitions 



t = h+t 2 



At = t 1 - t 2 



(33) 



In terms of these variables, some algebra yields the four 
possible decay amplitudes as 
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FF 

A n = ^-j= [(1 - 66 - V1V2) cos ±AAt 

+ (l + 66 + ^2)cos±AAAt 
-i(6 - 6) sin |AAt 
-i(6 + 6)sin§AAAf] e ~ 4At / 2 , 
A^ = -A /7 (6^-6,6^-6), 



./'./' 



= W -1 [(6V2 - 6^1) (cos i A At - cos ±AAAt) 



2V2" LVS1 " ~'^~"a* 
+*(Vi - V a ) sin ±AAt 
+i(^i + V 2 )sin±AAAt] e - jAt / 2 , 
Ayj = —Aff(F ^F,W^W-\b^ -6,6 -> -6), 

(34) 

where the dependence on £ and p\, p 2 is again under- 
stood. 

The expressions ([m]) are valid for CPT and T violation 
of arbitrary size and are independent of phase conven- 
tions. Nontrivial sensitivity to the sum and difference of 
6 and 6 is manifest. The corresponding decay probabili- 
ties are straightforward to obtain but are somewhat cum- 
bersome. They inherit the independence of phase conven- 
tions and the nontrivial sensitivity to 6 + 6- Since these 
factors depend on all four parameters Aa M for CPT vio- 
lation, appropriate analysis of experimental data for cor- 
related decays can provide four independent CPT tests. 

The type of analysis needed depends on the experi- 
mental situation. The remarks following Eq. about 
averaging and binning apply here, and there are also con- 
siderations specific to the case of correlated mesons. For 
example, if the quarkonium is produced at rest in the 
laboratory, perhaps by a symmetric collider, then the 3- 
momenta of the correlated mesons are equal in magnitude 
and opposite in direction. The sum 



6+6 = 2 7 (p)Aa /AA 



is then independent of Aa, so extracting an asymmetry 
sensitive to 6 + 6 yields a clean bound on Aao. Sim- 
ilarly, the difference 6 — 6 is independent of Aao, and 
binning in sidereal time permits bounds on the three com- 
ponents Aa. If instead the quarkonium is produced in an 
asymmetric collider, the two 3-momenta of the correlated 
mesons are not back-to-back in the laboratory frame, so 
6 ±6 are both sensitive to all components of Aa M . Four 
independent measurements of CPT violation can again 
be extracted. 

Many of the interesting features can be illustrated in 
the approximation of small 6> 6> f° r which the expres- 
sions simplify to some extent. This approximation is cer- 
tainly valid for the K system, and the recent results from 
OPAL, DELPHI, and BELLE Q imply it is also valid 
for the Bd system. The situation for the D and the B s 
systems is less clear, with large CPT violation remain- 



ing experimentally admissible, but many of the following 
considerations still apply. 

Consider for definiteness the double decay into //. To 
leading order in 6 and 6, the decay probability P^j is 



P f j = Pfj(t, At, i, pi,p 2 ) 

= j\FF\ 2 e~ 7 */ 2 {cosh±A 7 At + cos AmAt 
-Re (6 +6)sinhiA 7 Ai 
-Im (6 + 6) sin AmAt 
+2Im [(6 - 6) cos(±AA* At) sin(±AAt)] }. 

(36) 

This expression shows the combination 6 + 6 is associ- 
ated with an odd function in At, while 6 — 6 is associated 
with an even function in At. This distinction allows the 
separate extraction of 6 ±6- As an explicit example, 
the case of the sum 6 + 6 is treated here. 

In typical experimental situations for the correlated 
double-meson decay, the time sum t is unobservable but 
the difference At can be used as a fitting parameter. It is 
therefore appropriate to work with an integrated proba- 
bility T £j(At,t, pi, P2) obtained by integrating the prob- 
ability (pi) over t: 



T f j(At,i,p 1 ,p 2 ) 



A/ 



dt F /7 (t,At,t,p!,p 2 ). (37) 



An asymmetry ^4 CPT jj sensitive to the sum 6 + 6 °f 
parameters for CPT violation can then be defined as 



_ T f j(At,i,p 1 ,p 2 ) - Y f -j(-At,i,pi,p 2 ) 
r / j(At,t,pi,p 2 ) + T f j(-At,i,p 1 ,p 2 ) 

Calculation gives 



(38) 



(35) a 



L CPT,// - 

-Re (6 + 6) sinh \ A 7 At - Im (6 + 6) sin AmAt 
cosh ±A 7 At + cos AmAt 

(39) 

which is valid to lowest order in CPT-violating quanti- 
ties. For the Bd system, this expression generalizes the 
asymmetry obtained ||l5|| under the assumption of con- 
stant parameter for CPT violation and used to place the 
recent experimental limits on CPT violation at BELLE 
& 

For quarkonia produced in a symmetric collider the 
asymmetry ( |39| ) depends only on Aa because the sum 
6 + 6 is given by Eq. ( |35| ) . There is therefore no varia- 
tion with t, and the line spectrum of the mesons implies 
there is also no variation with p\ = —p 2 . In this case, a 
direct fit to the variation with At provides a bound on 
Aao- 
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In contrast, for quarkonia produced in an asymmetric 
collider the asymmetry ([59j) depends on all four parame- 
ters Aa M and also varies with t and pi, p 2 . For any given 
situation, forming an asymmetry of the type ( |38| ) after 
averaging Eq. ( |36| ) over suitable combinations of the vari- 
ables t, pi , p 2 permits the extraction of four independent 
CPT bounds, one for each parameter Aa M . Independent 
tests of this kind for the Bd system should be feasible at 
both BaBar and BELLE, where the quarkonia are pro- 
duced in asymmetric collisions and the meson pairs are 
boosted in the laboratory frame. 



VI. SUMMARY 

This work has studied some aspects of tests of CPT and 
Lorentz symmetry using neutral- meson oscillations. A 
formalism has been adopted for the treatment of arbitrar- 
ily large indirect CPT and T violation in the K, D, Bd, 
and B s systems that is phase-convention independent. It 
involves a real parameter w for T violation and a com- 
plex parameter £ for CPT violation. An expression for 
the latter, given as Eq. (|T^), is derived from the general 
Lorentz- and CPT-breaking standard-model extension. 
This equation reveals that CPT observables can vary 
with the magnitude and orientation of the meson mo- 
mentum and hence also with sidereal time. To illustrate 
some of the implications for experiment, transition am- 
plitudes, decay probabilities, and sample CPT-sensitivc 
asymmetries for semileptonic decays are derived. Both 
uncorrelated and correlated mesons are considered, and 
some consequences for experiments are described. 

The analysis shows that four independent experimental 
bounds are required to bound CPT violation completely 
in any single neutral-meson system. Since these param- 
eters may differ between systems, separate experimental 
analyses are required in each case. No bounds are avail- 
able in the D or B s systems as yet. Certain combinations 
of the four key parameters Aa^ have been constrained in 
the K and Bd systems by recent experiments |l^ , |l6| , but 
no definitive analysis has yet been performed. Obtaining 
a complete set of linearly independent measurements in 
any of the meson systems has the potential to offer our 
first glimpse of physics at the Planck scale and would in 
any case provide crucial experimental information on the 
existence of CPT and Lorentz violation in nature. 
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APPENDIX A: STANDARD FORMALISMS 

This appendix lists a few key properties of five stan- 
dard formalisms for indirect T and CPT violation. All 
these can be traced to early work several decades ago in 
the context of the K system || . For most of these stan- 
dard formalisms, several closely related variants exist in 
the literature, but for definiteness only one of each type 
is presented here. 

The MY formalism sets 

/ Mu - \%c n m 12 - §*r 12 \ 

A = M - \iT = . (Al) 

V M i*2 - ¥ V h M 22 - \iT 22 j 

The off-diagonal quantities are all phase-convention de- 
pendent. The parameter for CPT violation is the com- 
bination (Mn — M22) — — T 22 )/2. The parameter 
for T violation is \(Mf 2 - iT* 12 /2)/(M 12 - iT 12 /2)\. The 
masses and decay rates are given by 

A = (M 11 + M 22 )-iz(r u + r 22 ), 

AA = 2[(M 12 - \iT l2 ){M* 12 - \iT* l2 ) 

+ i [(Mn - M 22 ) - *i(r u - r 22 )] 2 ] 1/2 , (A2) 

where the definitions in Eq. (|J) are understood to hold. 
The DEiE 2 E 3 formalism sets 



A = 



-iD + E 3 E 1 - iE 2 



Ex + iE 2 -iD - E 3 



(A3) 



All off-diagonal quantities are phase- convent ion depen- 
dent. The parameter for CPT violation is E 3 . The 
parameter for T violation is i{E\E 2 — E\E 2 ). The 
masses and decay rates are given by A = —2iD, AA = 
2 VW 



E 2 



The DE9<j) formalism sets 



A = 



-iD + E cos 6 E sin ( 



EsmOe^ -iD -E cos ( 



(A4) 



The parameter <j) is phase-convention dependent. The 
parameter for CPT violation is cos#. The parameter for 
T violation is | exp(i(^>)|. The masses and decay rates are 
given by A = -2iD, AA = 2E. 

There are also formalisms that are introduced in terms 
of the relationship between the strong-interaction eigen- 
states P°, P° and the physical eigenstates P a , Pb- A 
general one is the pqrs formalism, which sets 



\P a )=p\P°)+q\n, 
\P b )=r\P°)- S \P0), 



(A5) 



where p, q, r, s are complex parameters. In this formal- 
ism, one can show 
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A 



( X(ps + qr) 
+AX(ps — qr) 



2{p± + qi) 



2AXpr \ 



\ 2A\qs 



X(ps + 
-AX(ps 



qr) 



Sanda, Eur. Phys. J. C 15, 1 (2000); I.I. Bigi, to appear 



(A6) 



The complex parameters p, q, r, s are all phase- 
convention dependent. They are also substantially re- 
dundant, since only three of their eight real components 
have physical meaning. The normalization conventions 
for the wave functions represent two degrees of freedom, 



often fixed by the choice \p\* 



= 1. The 



remaining three unobservable degrees of freedom are the 
absolute phases of |P a ) and \Pb) and the relative phase 
of \P°) and \P°). The parameter for CPT violation is 
(ps — qr). The parameter for T violation is \pr/qs\. The 
masses and decay rates are additional independent quan- 
tities, taken here as A, AA. 

The eS formalism |l5| is widely adopted for the K sys- 
tem. It can be regarded as a special case of the pqrs 
formalism. For arbitrary-size T and CPT violation, the 
eS formalism can be defined as 



\P a ) 

in} 



(i 



■S)\P°) 



a 



S)\po) 



(1 



■S)\P°) 



\e + 8\*) 



■S)\P°) 



x/2(l + |e-£| 2 ) 



(A7) 



In this formalism, A is given by Eq. (A6) with appropri- 
ate substitutions for the parameters p, q, r, s in terms 
of e, S, obtained from Eq. (A7). Both e and 5 depend 
on phase conventions. Nonzero values of e and S charac- 
terize T and CPT violation, respectively. For the special 
case of small e and 5, which is a good approximation in 
the K system, one can show 



A«i 



A + 2AA<5 AA(l + 2e) 
AA(1 - 2e) A - 2AX5 



(A8) 



Even within this approximation e is phase-convention de- 
pendent, although 5 is not. The parameter for T violation 
can then be taken to be Ree, for example. The masses 
and decay rates are independent quantities and here are 
specified by A, AA. 
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